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FINITE-THRUST ESCAPE FROM AND CAPTURE INTO 
CIRCULAR AND ELLIPTIC ORBITS 
by Edward A. Willis, Jr. 

Lewis Research Center 

SUMMARY 

Finite-thrust escape and capture trajectories are considered which lead from a cir- 
cular or elliptic orbit to a specified hyperbolic excess velocity vector located at the 
gravitational "sphere of influence. " Each trajectory consists of a powered maneuver 
initiated at the parking orbit and a coasting arc from the cutoff point to the sphere of in- 
fluence. The powered maneuver is accomplished with constant, continuous thrust and con- 
stant jet velocity. Tangential steering and optimum thrust- initiation points are used to 
obtain very nearly minimum characteristic velocity increments. 

These near-optimal finite-thrust trajectories are analyzed in dimensionless terms 
and compared with equivalent impulsive ones. The minimal characteristic velocity ratio 
(that is, the ratio of the actual propulsive effort, /a dt, to the equivalent impulsive veloc- 
ity increment), the optimum characteristic central angle (measured from a reference 
position on the parking orbit to the asymtotic direction at the sphere of influence), and 
the optimum initial true anomaly for departing from an elliptic orbit are found for a wide 
range of planetocentric trajectories. These results are presented as functions of the 
dimensionless hyperbolic velocity, acceleration, and jet velocity parameters. The di- 
mensionless arguments may be easily scaled to any particular case of interest - for ex- 
ample, a specific planet, parking orbit, propulsion system, and hyperbolic velocity. 

The data ranges were chosen primarily to correspond with high- to medium-thrust sys- 
tems, such as chemical rockets and various types of nuclear rockets. These results 
may be conveniently used in combination with ballistic interplanetary trajectory calcula- 
tions to obtain a realistic geometric description and accurate propellant fractions for a 
great variety of interesting and potentially important space flight missions. Examples 
illustrating the use of these data for actual mission problems are also presented. 



INTRODUCTION 


The dynamic and geometric effects of finite vehicle acceleration must be accounted 
for in order to derive accurate propellant fractions (and hence vehicle weight) from bal- 
listic interplanetary trajectory data (such as refs. 1 to 3). The velocity increments ob- 
tainable from impulsive calculations are useful in preliminary studies but do not always 
lead to accurate propellant fractions for vehicles with realistic thrust levels. This is 
because of the following: 

(1) A reduction in propulsive efficiency (also called a ’’gravity loss”) occurs when an 
impulse is replaced by a finite-thrust maneuver. As a result, the characteristic propul- 
sive velocity increment (AV ch ) and propellant fraction increase. 

(2) There is often a considerable difference between the geometric properties of im- 
pulsive and finite-thrust trajectories. This, in turn, may complicate the problem of 
matching heliocentric and planetocentric trajectory segments. 

The efficiency penalty due to finite thrust may be accounted for by applying a correc- 
tion factor (also termed characteristic velocity ratio) to the impulsive velocity increment. 
Accurate propellant fractions may then be obtained in a convenient and familiar way from 
the classical rocket equation. Prior discussions of this approach (e. g. , refs. 4 and 5) 
presented some typical velocity correction factors, but no geometric data, and consid- 
ered only circular parking orbits and escape maneuvers. An alternative approach 
(ref. 6) is to present the propellant fraction directly. In principle, this is equivalent to 
the correction factor approach but is less convenient to use because a greater amount of 
interpolation is required. A considerable amount of data of this nature is given in refer- 
ence 6, which is limited, however, to the case of circular orbits. The present report 
includes a wide range of velocity correction factors and geometric data for both escape 
and capture maneuvers and for elliptic as well as circular initial orbits. This report is 
intended both as a generalization of references 4 and 5 and to complement available and 
future ballistic trajectory data. The correction factors and geometric data presented 
herein may be used together with ballistic interplanetary data to obtain accurate propel- 
lant fractions and realistic overall geometry for many missions of practical interest. 

Trajectories that use a single burning period followed by coasting to the sphere of 
influence are analyzed herein. The powered maneuver begins at an optimum point on the 
initial orbit and proceeds with constant thrust and constant jet velocity until the hyper- 
bolic excess velocity attains a prescribed magnitude. After cutoff, the vehicle continues 
along a Keplerian arc to the sphere of influence. The final asymptotic direction of the 
hyperbolic excess velocity vector is left open. For such trajectories, it is shown that 
tangential steering yields very nearly the minimum characteristic velocities and propel- 
lant fractions. 

Characteristic velocity ratios, optimum initial true anomalies, and final asymptotic 
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directions are derived from a study of these unconstrained trajectories. The results are 
presented as functions of the initial orbit elements and the dimensionless hyperbolic ve- 
locity, jet velocity, and acceleration parameters. The significance of each dimensionless 
variable is discussed, and scaling rules are developed to apply the generalized data to 
any specific case of interest. As a sample application, the scaling procedures are used 
to generate a set of specific working data curves that describe typical nuclear rocket ma- 
neuvers near Venus, Earth, and Mars. Further examples illustrate the use of the pres- 
ent results in solving representative mission analysis problems such as the selection of 
an advantageous parking orbit radius and optimization of the initial thrust to weight ratio. 


SYMBOLS 


A acceleration 
a dimensionless acceleration 

e eccentricity of conic section 

F thrust 

f v characteristic velocity ratio or correction factor, AV^^/AV ^ ^ 

G acceleration due to gravity 
g dimensionless acceleration due to gravity 
H angular momentum 
h dimensionless angular momentum 
# auxiliary function, defined in equation (Al) 

I specific impulse, sec 

k structural mass fraction, see equation (43) 

M mass 
m mass fraction 
P semilatus rectum 
p dimensionless semilatus rectum 
R polar radius 

r radius ratio 

t time 
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u angle of attack 

V velocity 

v dimensionless velocity 

AV velocity increment 

Av dimensionless velocity increment 

a trajectory path angle 

7] efficiency 

9 trajectory central angle 

6 ch characteristic central angle 

/i planet gravitational constant 

v parking orbit true anomaly at thrust initiation (escape maneuver) or at thrust 

termination (capture) 

t dimensionless time 

^ adjoint variables, defined in equation (A2) 

()’ denotes differ entation with respect to r 

( ) denotes mean value 

Subscripts: 
act actual 

as acceleration sensitive 

bo burnout 

c circular 

cap capture 

ch characteristic 

cst coasting 

esc escape 

g gross 

i initial 

imp impulsive 

j jet 
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k, n general numerical indices 

SL payload 

m mean 

max maximum 

min minimum 

opt optimum 

p propellant 

pi planet 

po parking orbit 

pr propulsive 

ps propellant sensitive 

pwr powered 

ref reference 

sc spacecraft 

ts thrust sensitive 

°° sphere of influence 

® Earth 

$ Venus 

<f Mars 


ANALYSIS 

The initial mass of a space vehicle is a useful criterion for interplanetary mis- 
sion studies. A major step in computing Mj is to determine the propellant fraction m^ 
for each propulsive maneuver. The propellant fractions are derived in turn from a study 
of the flight trajectory, which is illustrated in typical form in figure 1(a). The flight 
path consists of alternate planetocentric arcs (view A) and heliocentric arcs (view B), 
which are matched at the spheres of influence as indicated in view A; that is, 

V °o, planetocentric trajectory = ( V sc " V pl) he ii 0 centric trajectory W 
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Either type of arc may contain a propulsive maneuver. For all but very low-thrust 
systems, the propulsive effort or AV^ of a heliocentric maneuver (used as a midcourse 
correction or, as in ref. 3, to reduce the total AV) may be determined accurately by im- 
pulsive methods. Values of m are then given by the familiar formula 

ir 

(2> 

Planetocentric trajectories (illustrated in figs. 1(b) to 1(d)) join a parking orbit and 
must produce a ’’hyperbolic excess velocity” vector V 00 (as defined by eq. (1)) at the 
sphere of influence. Relative to the parking orbit, is represented by its magnitude 
and direction 0 These two parameters together with the parking orbit elements de- 
fine boundary conditions for each end of the planetocentric trajectory. Impulsive maneu- 
vers to satisfy these boundary conditions (fig. 1(b)) can be easily computed, but they often 
lead to erroneous values of m . The effects of finite thrust and other neglected items can 
be lumped into a "correction factor" f v , which will be so defined that accurate propellant 
fractions for finite-thrust maneuvers (figs. 1(c) and 1(d)) can still be obtained from an ex- 
pression of the same simple form as equation (2); that is, 

/f AV im \ 

m P = 1 ■ (3) 

The quantity f y AV^ m p may be recognized as the actual propulsive effort or character- 
istic velocity AV^ for the finite-thrust maneuver. Thus f v is defined, for compu- 
tational purposes, as 

t = a, v= 1 . r°° - Thr Hg dt (4) 

V av imp aV im P y 0 Mass 

and is termed the characteristic velocity ratio. 

It is convenient to define AV. „ as the minimum AV which will produce the same 

imp 

magnitude of from the same parking orbit as does the finite-thrust maneuver (with- 
out regard to the value of 0^). As is well known, the control policy that minimizes 
AVjmp is to apply the impulse tangent to the parking orbit and at periapse if the orbit 
is elliptic. With this definition, AVj m p can be easily computed in closed form. Since 
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AV. is minimal, f ^ 1 (it will be shown later that f is also bounded below 3. 0 in 

imp v v 

the case of circular orbits). 

On the other hand, AV^ is obtained by numerical integration, thus accounting for 
finite thrust and jet velocity, possible geometric constraints on 0 and actual (as op- 
posed to optimal) control policies in addition to and the parking orbit elements. 
Consequently, f^ reflects these same factors. 

For impulsive trajectories, the characteristic angle consists only of a coasting arc 
0 imp as fi& ure 1(b) illustrates. When finite thrust is used, 0 ch consists of a powered 
central angle 0 pwr plus a coasting arc 0 cgt as indicated in figure 1(c). The sum of 
these two contributions is generally larger than the coasting angle 0 imp , which corre- 
sponds to an impulsive trajectory. Thus, the power-on point must be advanced in the 
finite-thrust case in order to attain a prescribed direction of in inertial space. 

The quantities f and 0^ define the essential properties of finite- thrust trajec- 
tories that begin or end in circular or elliptic orbits. The values of m^ obtainable 
from precomputed f v data by means of equation (3) are more accurate than those which 
could be obtained if m^ were plotted directly (as in ref. 6). This is because AV irn p 
is always maintained as the first approximation to AV ^ For the same reason, the 
present approach is often more convenient to use; in many cases, constant "typical** 
values of f^ can be chosen that will yield acceptably accurate values of without 
laborious interpolation. 


Assumptions 

The following assumptions have been introduced in order to simplify the numerical 
calculations: 

(1) As previously mentioned, an overall trajectory consists of successive two-body 
arcs joined at the planet’s gravitational sphere of influence. 

(2) Each arc lies in a plane determined by and the center of force, and the 
gravitational attraction along each arc is represented by an inverse square central force 
field. 

(3) The propulsion system operates with constant continuous thrust and constant jet 
velocity. The thrust is considered finite, but large enough that the necessary magnitudes 
of can be attained within the sphere of influence. The radius of the sphere of in- 
fluence is assumed to be much larger than that of the initial planetocentric orbit. 

(4) The sphere of influence boundary condition is taken to be the magnitude | V M | 
only. The characteristic angle 0 ch , which defines the orientation of the initial orbit rel- 
ative to V OT , is not constrained but is treated as a dependent variable. Therefore, for 
given V m and parking orbit elements, the factors that primarily affect 0 ch , namely the 
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power-on or -off point v (cf fig. 1(d)) and the steering control policy u(t) (i. e. , the 

po 

angle of attack program), can be chosen solely to minimize AV^ and hence f y . Tra- 
jectories of this kind will be termed ’’optimum-angle” trajectories herein. Such trajec- 
tories can always be used with a circular orbit since any necessary direction of V M in 
inertial space can be attained, without changing f , by selecting the appropriate power- 
on or -off point. The optimum angle can also be used for one maneuver based on a given 
elliptic orbit. If an elliptic orbit is used at the destination planet of a round-trip mission, 
it must be reoriented between arrival and departure in order to use optimum-angle tra- 
jectories for both capture and escape. Several effective methods of reorienting an elliptic 
orbit are described in references 7 and 8; these methods would permit optimum-angle 
trajectories to be used for both escape and capture maneuvers. 

(5) Tangential steering (i. e. , zero angle of attack) is used to approximate the optimal 
steering policy. If 0^ is not constrained, tangential steering yields the minimum pro- 
pulsive effort in the limit of impulsive thrust. In reference 9 it was shown that, with 

0 . not constrained, the penalty for using tangential rather than variational steering for 
cn / o \ 

constant acceleration escapes to parabolic energy = 0) from a circular orbit is quite 

small (under 1 percent) regardless of the mangitude of the acceleration. In appendix A, 

the analysis of reference 9 is generalized to include constant thrust, constant jet velocity, 

hyperbolic burnout energy (V^ > 0), and elliptic as well as circular orbits. It is shown 

that the maximum penalty for tangential steering is still under 1 percent of AV^ under 

these more general conditions, and it decreases rapidly as elliptic orbits of increasing 

eccentricity are considered. 


Equations of Motion 

When the foregoing assumptions are used, a simple set of differential equations may 
be constructed to describe the rocket's motion. 

Coordinate system an d dimensionless variables. - The trajectory is represented in 
plane polar coordinates as indicated in figures 1(b) to 1(e). The position of the rocket R 
is described by the polar radius R and central angle 0 (fig. 1(e)); the velocity V is 
denoted by its magnitude V and path angle a relative to the local horizontal. Similarly, 
the acceleration A is described by its magnitude F/M and angle of attack u (relative 
to V). Gravity G is represented by its magnitude p/R 2 and is directed toward the 
center of force. 

In order to obtain general results, it is necessary to normalize the variables just 
described. The resulting dimensionless variables (denoted by lower case symbols) are 
obtained by dividing radii, velocities, accelerations, and time, respectively, by the 
radius, circular velocity, local gravity, and circular radian period corresponding to a 
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reference position on the initial orbit. Polar angles (0) are also measured relative to the 
reference position. For circular orbits, the reference position is the power-on or -off 
point, but for elliptic orbits the periapse is the reference point. The dimensionless vari- 
ables, which will be used henceforth, are related to their dimensional counterparts by 


rfr) = - R &- 


R 


ref 


V(T) 

V c, ref 


a( T ) = 

f//R: 


ref 


T= ^c 1 ref t 

R ref 

G 1 


g = 


M/ R ref ^ 


V 


G © 1 


c, ref 


v 2 = 

CO 


c ? ref/ 


(5a) 

(5b) 

(5c) 

(5d) 

(5e) 

(5f) 

(5g) 


Basic equations . - With the normalization just described, the first-order vector 
equations of motion are 


t 

r 


= v 


V = a +g 


The components of equation (6) for r, v, a, and 6 are then 


(6a) 

(6b) 
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r* = v sin a 


(7a) 


v* = a cos u - 


sin a 


(7b) 


a t = a sinu + / 2 _ 1\ cos a 
v \ r/ rv 


= — cos a 
r 


(7c) 

(7d) 


The characteristic velocity AV^ which is needed to define f v is given by 

AV^ h = a(r) 


( 8 ) 


Since the vehicle operates with constant thrust and jet velocity, the acceleration a(T) is 
given, for an escape maneuver, by 


a( T ) 


po 


esc 1 / 

1 - a P 0 T/v i 
= 0 


0 < t < r 


bo 


< r (coasting flight) 


(9a) 


where ap Q = a^, the initial acceleration at the power-on point. 

Capture maneuvers, on the other hand, are integrated backward in time, that is, be- 
ginning at the desired power-off point on the orbit. This is done only for computational 
convenience (to avoid the iterative search that would otherwise be required) but has the 
consequence that the acceleration term for capture maneuvers is given by 


a(T) 


= 2^ 0 < T < T. 

cap 1 + V T/v i 


= 0 < t (coasting flight) (9b) 

In this case ap Q is to be interpreted as the acceleration at burnout rather than the initial 
one. For either escape or capture maneuvers, the initial and burnout accelerations are 
related by 
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(10a) 


a bo " 


a i 


1 - m_ 


1 - 


Vbo 


or 


a bo 

a. = — 

1 , a bo T bo 


“bo 



M- 


bo 


(10b) 


Boundary conditions . - In terms of the present dimensionless variables, the initial 
(r = 0) conditions for a circular orbit are 


r(0) = 1 

(Ha) 

v(0) = 1 

(lib) 

a(0) = 0 

(He) 

0(0) = 0 

(lid) 


For an elliptic parking orbit, the initial conditions are given in terms of the eccentricity 
of the orbit ep Q and the true anomaly of the power-on or -off point v^ Q by the well- 
known conic equations 
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1 + e__ 

r(0) = - 1)0 


1 + e po cos G(°) 


v(0) = 


L r (o) 


- (» - V 


1 1/2 


Q'(O) = tan 


-1 


| e po r (°) sin 0(°)| 


1 + e. 


po J 


0 ( 0 ) = v 


po 


(12a) 


(12b) 


(12c) 


(12d) 
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With these initial conditions and a known form for the steering control law (for ex- 
ample, u(t) s 0), the equation of motion may be integrated (numerically) to obtain r, v, 
a, and 6 as functions of r with v , e , v . a , and v- as parameters. Thrust is 

pu pu 00 pu J 

terminated at the burnout time r^, at which time the desired energy or hyperbolic veloc- 
ity has been attained (cf eq. (1)), that is, when 


2 / \ 2 , \ 2 2 
v J T bo> = v < T bo> ' -7—, = v co 


bo' 


’< T bo> 


= |Vsc ~ V P lL (13) 


at sphere of influence 


V 2 

v c, ref 


Once t^ o has been determined, equations (8) and (9) may be evaluated as a definite inte- 
gral; that is, 


Av ch( T bo> A 


C T bo 

t/n 


a(r)dr = Vj In 



(14) 


The properties of the remainder of the trajectory, that is, the coast from burnout to the 
sphere of influence, may finally be calculated in closed form from the well-known conic 
equations. 


Trajectory Solutions 

There are no known general closed- form solutions to the equation of motion described 
above. As has been previously mentioned, accurate values of Av^, 0 and other tra- 
jectory parameters can only be obtained by numerical integration. Impulsive solutions, 

which can be obtained in closed form, are required to define values of Av. and thus 

imp 

of fy. Low-thrust limiting solutions are also studied to derive an upper bound for f y 
and to obtain insight into the behavior and significance of this function. 

Impulsive thrust . - For tangential steering (u(r) = 0), only the velocity equation (7b) 
contains a. For sufficiently high thrust, and/or for very low local gravity, the dimen- 
sionless acceleration a tends toward infinity during the powered maneuver, but the dura- 
tion of the maneuver tends toward zero. In the impulsive limit, v changes discon- 

tinuously across the maneuver, but r, a, and 0 remain constant. The velocity equa- 
tion (7b) can then be integrated in closed form to yield the "classical rocket equa- 
tion" (2). The appropriate value of Av imp (= v bQ - v^ may then be derived from equa- 
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tions (12b) and (13) when it is recalled that the impulse is not only tangential but is to be 
applied at the periapse if the orbit is elliptic; that is, 


Av 


-(4 + i) 1/, -(i + 0 


imp 


1/2 


pOv 


(15) 


This is the value which enters the definition of f y (eq. (4)). 

The remainder of the planetocentric trajectory consists of a coasting arc (i. e. , a hy- 
perbola) whose periapse is at the impulse point and whose asymtotic branch extends to 
the sphere of influence (recall fig. 1(b)). Since the sphere of influence is considered very 
large compared to the parking orbit, the characteristic angle 0 ch is well approximated 
by the asymtotic true anomaly of the hyperbola; that is, 


0 • „ — cos 

ch, imp 


-1 



(16) 


where 


e bo" 


+ 1 


(17) 


Finite thrust. - The powered portion of the trajectory is integrated numerically to 
obtain values of r, v, a, and 0 at the burnout time r^ 0 , in addition to a value of AV^ 
for use in equation (4). The remainder of the trajectory, that is, the coast from the burn- 
out point to the sphere of influence, is again a conic arc whose elements may be computed 
in terms of the burnout conditions r^, v^ Q , a^, and The semilatus rectum p^ 

is given by 

Pbo ■ r bo v bo cos2 “bo < 18a) 


and the eccentricity e bQ by 


e bo 



(18b) 


Then the coasting contribution (0 j.) may be computed from the conic equation and added 

to 0. _ to obtain 
bo 


14 



e ch - 6 bo + cos 


-1 


J_\ - cos" 1 ' 1 C>bo/lbo 


(19) 


e b Q y 


bo 


Very low thrust , - A useful approximate solution may be obtained in closed form for 
maneuvers with very low but finite thrust. It is assumed that 


a (t) « — (r) 

2 

r 


< 0 


(20a) 


but that 


a(T) » — (r) 
r 2 


>« 


(20b) 


That is, escape (v^ = o) is attained at a very great radius, and further energy is then 
added in a field-free region. If the initial orbit is circular it will remain nearly circular, 
that is, a spiral of low pitch. This means that a is small and moreover that a' » 0. 
Equation (7c) then yields, for tangential steering, the condition 


2 „ 1 
y « — 

r 


( 21 ) 


Combining this relation with equations (7a) and (7b) yields 


( Av ch)' = v ’ + ' 


.3/2 


(22 a) 


or 


d(Av . ) = dv + d ( — 

wy 


(22b) 


Integrating term by term and applying initial condition (11) and the terminal conditions 


r ( T bo ) - 00 
v < T bo) ^ v - 


(23) 


yield the result that 
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Av ch, very low thrust, ~ v °° 
circular orbit 



+ 1 


(24) 


so that 


very low thrust, 
circular orbit 


V CO + 1 


/ 2 \l/2 

WL + 2 ) - 1 


( 25 ) 


It will be shown in the next section that Av ch is directly proportional to the change 

in v 2 that occurs in going from the initial orbit to the burnout point; that is, 

2 


Av C h a 


(v» bo ~ v °° po)‘ su S§ es<; s that Av^ for an elliptic orbit departure should 

be approximately the same as for a departure from a circular orbit whose energy is equal 
to the energy of the ellipse. Since 


v 2 - _ i _ 2 

v °°, elliptic po “ e po “ v =°, equivalent 


-1 


(26) 


circular po c ’ ^valent 


the radius of this equal- energy circular orbit is 


P — _ 

c, equivalent j _ e 

po 

If (27) is used in place of (11a), equation (24) becomes 


(27) 


AV ch, ellipse, =Vbo +(1 -V l/2 (28) 

low thrust 

and the corresponding value of f y is 


v„ i — + (1 " O 


°°, bo 


1/2 


po' 


*v,very low thrust, ~ j/2 

elliptic po (v 2 ? bQ + 2) - (1 + e pQ ) 1/2 


(29) 


It may be seen by inspection of equation (25) that f y for a circular orbit has a maxi- 
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mum value of 3 where v^ = 0. 5. This is the desired upper bound. For elliptic orbits, 
it may be seen from (29) that the upper bound is an increasing function of e po . If 
e^ = 0. 9, for instance, the maximum value of f y is about 10. 6 (rather than 3) and oc- 
curs when v^ « 0. 15. It will be demonstrated in the RESULTS that equations (25) and 
(29) not only yield valid upper bounds for f^ but also provide an accurate qualitative de- 
scription of how f y varies with v^ for low-thrust systems. 

M ean value solutions . - A formal solution for Av^ useful in interpreting the behav- 
ior and significance of f y may be derived by converting the velocity equation (7b) into an 
energy equation. Multiplying by v and using equation (7 a) yield 

vv' + — = va cos u (30a) 


or (recalling eq. (13)), 


1 

2 



* = a(v cos u) 


(30b) 


This may be formally integrated, using the law of the mean, to obtain the symbolic result 


2 

oo. p 0 


Ay. g "■»» ^ 

c 2(vcosu) 


( 31 ) 


where (v cos u) is the mean tangential velocity obtained by averaging v cos u with re- 
spect to Av^. When the same argument is applied to an impulsive trajectory, it is seen 
that 


, ± Av eh _ ( v cos u ) imp 

_av imp < vcosu > actual 


(32a) 


That is, f may be expressed as the ratio of the mean tangential velocity for an impulsive 
maneuver to that of the actual finite-thrust maneuver. In principle, f^ may be deter- 
mined by comparing the velocity and steering histories of actual and impulsive maneuvers. 

o 

(Examples of this will be given in the DISCUSSION. ) For a given value of Av ot , the 
parameters a pQ , Vj, and v only affect by their influence on v. 

On the other hand, the steering control u(t) enters f y explicitly, as is clearly 
shown in (31) or (32a). It is evident that tangential steering (u(t) = 0) will produce the 
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maximum instantaneous value of v cos u. Improvements in the mean value ( v cos u) by 
means of variational steering can only come about by increasing v at the expense of 
cos u. These opposing trends explain the fact, demonstrated in appendix A, that there is 
little to be gained by using variational rather than tangential steering for optimum-angle 
trajectories. If tangential steering is used, (32a) becomes 

< v) . „ 

fy a 1522- (32b) 

actual 

Thus, the selection of the power-on point v to minimize f^ is in accordance with the 
intuitive principle that energy can be added most efficiently in a region of high velocity. 

In fact, it is easy to see that ^ is a direct measure of the mean propulsive efficiency 
defined by 


V 


pr,m 


\ M P V J 2 


Energy added to final mass 


Useful energy liberated by propulsion system 


(33) 


Then, the ratio of 77 pr m for an actual maneuver to that for the equivalent impulsive 
manuever is (assuming equal M^) 


^P r ’ m ^ actual . M p, imp , 1 
Ov, m) imp M p, actual f P 


If equations (2) and (3) are used for M , 1/f may be expressed as 

it r 

1 Av lm,/ v l> - 1 

f p ex P( iv lm p/ v j) - 1 


fj , V - 1 Mmp) 

2f - 3 Av. \ 
i . v / imp) 

— 

> 

2 u > 

L 6 UJ 


> 

- 

J 


( 35 ) 


Therefore, 
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( 36 ) 


1 _ (V , my imp 
*p ^pr, m) actual 

to the first order in (^ v j m p/ v j)- Thus, f y varies inversely with the mean propulsive ef- 
ficiency. 

f y as a "gravity loss". - It is evident from equation (31) that f contains both dy- 
namic effects (i. e. , Av^) and geometric effects «v cos u)), neither of which depends ex- 
plicitly on the presence or absence of a gravitational field. In the general case, f can 
represent "steering losses" as much as gravity losses. For the optimum-angle trajec- 
tories which are of primary interest, however, there are no steering losses as such; 
therefore, f may be interpreted as a gravity loss in the sense that it accounts for dy- 
namic phenomena (reduction of ( v) and rj ) that are caused in turn by the gravity 

pr , 111 

field. As here defined and used, however, f y cannot be directly related to a kinematic 
gravity loss of the form 


Av 


gravity loss 


= f h ° sin SL 

Jo 1,2 


dr 


(37) 


When tangential steering is assumed in (7b), it may be seen that 


Av gravity loss Av ch " Av 


(38) 


while the loss associated with f is 

v 


Av loss, f = Av imp( f v “ = Av ch - Av : 


imp 


(39) 


A definition of f based on equation (37) would not necessarily represent a loss of 
energy since it would apply to coasting (constant energy) orbits as well as to powered 
maneuvers. Furthermore, such an interpretation cannot be used to explain the fact that, 
in the low-thrust limit, f continues to increase beyond v^ = 0. At this point (r — °°) all 
of the gravity field has been left behind, but as may be seen from equation (25), f in- 
creases from 2. 414 to 3. 000 as v^ increases from 0 to 0. 5. It will be shown in the 
DISCUSSION that this typical behavior can be adequately explained in terms of (v) . 
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RESULTS 


The equations of motion and auxiliary formulae indicated in the preceding section 
were programed for numerical integration on a high-speed digital computer. Values of 
fy and 0 ch were then obtained which correspond to a wide range of propulsive maneuvers. 


Dimensionless Data 


The quantities f y and 0^, which describe the overall dynamic and geometric char- 
acteristics of optimum -geometry trajectories, are plotted against the dimensionless hy- 

2 

perbolic velocity parameter v w . Both escape and capture maneuvers are considered 
over a range of values of the dimensionless variables and Vj. 

Results pertaining to maneuvers leading from or to a circular orbit are presented in 
figure 2. The dependent variable (i. e. , f y or 0 ch ) and the range of dimensionless agru- 
ments covered by the figure are listed in table I. It has been previously noted that any 
desired orientation of v M inertial space can be attained without changing f v simply by 
selecting the appropriate thrust- initiation point on the ciruclar orbit. Consequently, the 
value of is of little concern for circular orbits and 0 ^ data are presented for only 
two cases - figures 2(b) and 2(i). 

As indicated in table IE (p. 29), the corresponding results for elliptic orbits (with 
e^ = 0. 9) are presented in figure 3. In this case, 0^ defines the required orientation of 
the ellipse with respect to v M . It may be recalled by again referring to figure 1(d) that 
reference conditions for an elliptic orbit are defined by the periapse radius vector. The 
characteristic central angle 0 ^ is measured from periapse regardless of the initial 
power-on point (u ). The reference impulsive Av (eq. (3)) used in defining f is com- 

po V 

puted at periapse even if the actual trajectory begins at some other position (v * 0) on 

po 

the ellipse. 


Conversion of Dimensionless Data to Dimensional Form 

The preceding dimensionless results can be applied in any case of physical interest. 
That is, the reference values R re f, V c etc. are known once the planet and the park- 
ing orbit have been selected, while and/or AVj m p are known as a result of inter- 
planetary calculations. The dimensionless variables v^,, ap Q , and Vj corresponding to 
the maneuver may then be found by means of equation (5) and used with figures 2 and 3 to 
determine values of f and 0^. 

In some interplanetary studies (ref. 1), the results are presented in the form of 
values of v M for each trip. Other studies, such as references 2 and 3, present the 
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TABLE I. - KEY TO POWERED ESCAPE AND CAPTURE TRAJECTORY DATA, 
CIRCULAR PARKING ORBITS, AND DIMENSIONLESS ARGUMENTS 


Figure 

Depen- 

Maneuver 

Jet 

Dimensionless 

Hyperbolic 


dent 


velocity 

acceleration 

velocity 


variable 


parameter, 

at parking 

parameter, 




v j 

orbit, 

v i 




a po 


2(a) 

f v 

Escape and capture 

oo 

0-5.0 

-1 - 20 

2(b) 

d ch 

Escape and capture 

OO 

0-5.0 

-1 - 20 

2(c) 

f v 

Escape 

5 

0-5.0 

-1 - 20 

2(d) 

f v 

Capture 

5 

0-5.0 

-1 - 20 

2(e) 

f v 

Escape 

2. 5 

0-5.0 

i 

w 

i 

to 

o 

2(f) 

f v 

Capture 

2. 5 

0-5.0 

-1 - 20 

2(g) 

f v 

Escape 

1.0 

0-1.0 

-1 - 10 

2(h) 

f v 

Capture 

1.0 

0-5.0 

-1 - 10 

2(1) 

0 ch 

Escape and capture 

1.0 

0.01 - 1.0 

-1 - 10 

2(j) 

f v 

Escape 

0. 5 

0 - 0. 75 

-1 - 10 

2(k) 

f v 

Capture 

. 5 

0-5.0 

-1 - 10 

2(£) 

f v 

Escape 

0.25 

0 - 0. 4 

-1 - 10 

2(m) 

f v 

Capture 

.25 

0-10.0 

-1 - 10 



Characteristic velocity ratio. 



~ ~ O t u xu x*t 

Hyperbolic velocity parameter, v£ 

(a) Velocity correction factor; jet velocity parameter, ®. (b) Trajectory deflection angle,- jet velocity parameter, «. 

Figure 2. - Properties of escape and capture maneuvers. Circular parking orbit 
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Characteristic velocity ratio, f v ■ AV rh /AV 







Hyperbolic velc 

(g) Velocity correction factor; escape maneuver; jet velocity param- 
eter, 1.0. 


Figure 2. 







Trajectory deflection angle, 0 C | 1 , deg 


Escape 

Capture 



Dimensionless a. 1,0 f 

-acceleration at" -- 
_ parking orbit, j < 1.4 

a I 

P° ^ 






O 

+-• 

03 



>» 


*..05 

> 

1.1 


U 

1.08 


*t/5 


, . 1 

CD 

■#-» 

1.06 


O 

£ 


. r .L 

CO 

sz 

CJ> 

1.04 




“Dimensionless 
-acceleration at 
— parking orbit, 

~ 3nn 



Hyperbolic velocity parameter, v£ 


(i) Trajectory deflection angle; jet velocity parameter, 1.0. 


(j) Escape maneuver; jet velocity parameter, 0. 5 


Figure 2. - Continued. 









TABLE II. - KEY TO POWERED ESCAPE AND CAPTURE TRAJECTORY DATA 


[Elliptic parking orbits; parking orbit eccentricity, 0. 9; optimum values of op- 
timum orbital true anomaly, y po ; optimum trajectory deflection angle, 0 ch *] 


Figure 

Dependent 

variable 

Maneuver 

Jet 

velocity 

parameter, 

v i 

Dimensionless 
acceleration 
at parking 
orbit, 

a po 

Hyperbolic 

velocity 

parameter, 

3(a) 

f v 

Escape and capture 

OO 

0-1.0 

0-20 

3(b) 

0 ch- V 

Escape and capture 

OO 

0-1.0 

0-20 

3(c) 

f v 

Escape 

5.0 

0.01 - 1.0 

0 - 20 

3(d) 

e ch’ p o 

Escape 

5.0 

.01 - 1.0 

0-20 

3(e) 

f v 

Capture 

5.0 

.01 - 1.0 

0-20 

3(f) 

0 ch’ v po 

Capture 

5.0 

.01 - 1.0 

0 - 20 

3(g) 

f v 

Escape 

2.5 

0.01 - 1.0 

0-20 

3(h) 

0 ch>Po 

Escape 

2.5 

.01 - 1.0 

0-20 

3(i) 

f v 

Capture 

2. 5 

.01 - 1.0 

0-20 

30) 

^ch’ ^po 

Capture 

2. 5 

.01 - 1.0 

0-20 

3(k) 

f v 

Escape 

1.0 

0.01 - 1.0 

0-12 

3(1) 

^ch’ p o 

Escape 

1.0 

.01 - 1.0 

0-12 

3(m) 

f v 

Capture 

1.0 

.01 - 1.0 

0-12 

3(n) 

^ch’ ^po 

Capture 

1.0 

.01 - 1.0 

0-12 

3(o) 

f 

V 

Escape 

0. 5 

0.01 - 1.0 

0-12 

3(p) 

^ch’ ^po 

Escape 

. 5 

.01 - 1.0 

0-12 

3(q) 

f v 

Capture 

. 5 

.01 - 1.0 

0-12 

3(r) 

e ch’ > 

Capture 

. 5 

.01 - 1.0 

0-12 

3(s) 

f v 

Escape 

0.25 

0.01- 0.2 

0-8 

3(t) 

e ch> > 

Escape 

.25 

.01 - .2 

0-8 

3(u) 

f v 

Capture 

.25 

.01 - 10.0 

0-8 

3(v) 

9 ch> V 

Capture 

.25 

.01 - 10.0 

0-8 



Characteristic velocity ratio, f v = AV 


co 

o 



Figured - Properties of escape and capture maneuvers. Elliptic parking orbit; parking orbit eccentricity, 0.9. 
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Hyperbolic velocity parameter, v£ 


(l) Optimum deflection angle and true anomaly, escape; jet velocity pa ram 
eter, 1.0. 


Figured - Continued. 
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Figure 3. - Continued. 









. equivalent impulsive velocity increment necessary to leave or enter an assumed 

circular parking orbit at 1. 1 planet radii. In figure 4(a), v^ is plotted as a function of 
AVim (* n miles /sec) for each planet of the solar system, with Rp Q = 1. 1 Rp^. Figure 
4(b) gives the conversion between AVj m p and v^ in EMOS (i. e. , in units of the Earth’s 
mean orbital speed, 18. 5058 miles/sec). In figure 4(c), a p0 is shown in terms of A po 
(in units of Earth surface gravity, 32. 17 ft/sec 2 ), again assuming = 1.1 R pl . The 
dimensionless jet velocity parameter Vj is shown in figure 4(d) as a function of specific 
impulse I (in sec). 

As an example, suppose that it is desired to attain a value of V M = 0. 2 EMOS from 
a circular Earth-centered parking orbit (R^ = 1.1 R ffi ) using a chemical rocket with 

I = 385 seconds and (= A^) = 9 feet 
per second 2 or 0. 27 G 0 g . From fig- 
ures 4(a) and (b) it may be seen that 
the corresponding value of AVj mp is 

2 

about 3 miles per second and v m is 
approximately 0.7. Also, a^. « 0. 3 

pu 

and Vj ~ 0. 5, as noted in figures 4(c) 

and (d). If these values (v^ = 0. 7, 

a = 0. 3, v- = 0. 5) of the dimension- 
po ] 

less arguments are used, the corre- 
sponding value of fy may be found in 
figure 2(j) (p. 25) - that is, f^ = 1. 037. 
Thus, the actual characteristic veloc- 
ity increment is AV^ « 1 . 037 x 3 or 
3. 11 miles per second, and the propel- 
(a) Hyperbolic velocity parameter. lant fraction is 




0 2. 5 5. 0 7. 5 10. 0 12. 5 15. 0 

Impulsive velocity increment, AV|, miles/sec 
(b) Hyperbolic excess velocity. 


Figure 4. - Scaling of dimensionless parameters for planets of the 
solar system. 


m p = 1 - exp 


-3. 11 x 5280 
385 x 32.17 


or 0.735. Exactly the same procedure 
would be used for a capture maneuver 
with a specified value of A^ Q (rather 
than Aj); when figure 2(k) is used, the 
value of fy for the arguments enumer- 
ated above is seen to be f y = 1. 230. 

On the other hand, if it is desirable 
to study capture maneuvers on the 
basis of given values of Aj, or to 
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Jet velocity parameter, V; Dimensionless acceleration, 



Thrust to weight ratio at orbit, Gpo 
(c) Acceleration parameter. 



Figure 4. - Concluded. 


study escape maneuvers with given 
values of A bo , it is necessary to solve 
equation (10) by iteration. For instance, 

to find f for the previous capture ma- 

* 2 
neuver with = 9 feet per second re- 
quires the following steps: 

(1) The v M and v^ are as given 
previously. 

(2) A trial value of ap Q (= a^) is 
assumed. 

(3) The f y is found in the data and 
a trial value of m^ is computed by 
using equations (3) and (15). 

(4) A new value of a (= a^) is 
computed from equation (10a) - compare 
with step (2). 

(5) Steps (2) to (4) are repeated using 
the new value of a^ from (4) until con- 
vergence is obtained. 

A similar procedure is used to find the 
characteristics of escape maneuvers 
with given values of A^ Q . 

Dimensional Data 


For studies which involve repetitive 
but similar calculations, it is convenient 
to perform the previous conversion, interpolation, and iteration in advance. This yields 
working charts that apply to a specific planet and parking orbit, and to ranges of A^ 

(or A bo^ and I which are of particular interest. As an example, the preceding dimen- 
sionless circular orbit data have been used to generate a series of working charts that 
describe nuclear rocket maneuvers near Earth, Venus, and Mars (figs. 5, 6, and 7, re- 
spectively). Values of f^ are plotted against AV^ m p, assuming a circular orbit radius 
of 1. 1 planet radii for values of A^ ranging from 0. 1 to 1.0 G 0 and specific impulses 
of 900 seconds (solid curves) and 700 seconds (dashed curves). Part (a) of each of these 
figures represents escape maneuvers while captures are shown in part (b). It is noted 
that the captures are correlated on the basis of the initial acceleration Aj. 


41 




Characteristic velocity ratio, 


N> 



(a) Earth escape maneuver. Initial circular orbit. (b) Earth capture maneuver. Final circular orbit. 

Figure 5, - Characteristic velocity ratio at Earth. Constant tangential thrust; altitude, 1. 1 planet radii. 



Characteristic velocity ratio, f v = AV ch /AVj; 



(a) Venus escape maneuver. Initial circular orbit. (b) Venus capture maneuver. Final circular orbit. 

Figure 6. - Characteristic velocity ratio at Venus. Constant tangential thrust; altitude, 1. 1 planet radii. 
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Characteristic velocity ratio, f v * AV ch /AV| 


4 ^ 

4 *> 



Figure 7. - Characteristic velocity ratio at Mars. Constant tangential thrust; altitude, 1. 1 planet radii. 



DISCUSSION 


The preceding results define the essential dynamic and geometric properties of a 

very wide range of optimum-angle finite-thrust trajectories. This section is intended to 

interpret the behavior of these results by examining the role of each of the independent 

2 

variables. The effects of v M , a^, and Vj are first considered in terms of circular 
orbits. The additional variables (v w and e^ associated with elliptic orbits are then 
discussed, and finally the circular and elliptic orbit results are compared. 


Circular Orbits 

o 

Effects of hyperbolic velocity parameter v^. - The characteristic velocity ratio f y 

is shown as a function of v^ (and a^) in figure 2(a) (p. 21) for escape or capture maneu- 
vers with Vj = °° (for instance) . The uppermost curve on this graph is the low- thrust 
limiting solution expressed in equation (25). This curve gives the upper bound for ^ as 
a function of v^; the peak value of ^ = 3 occurs at v co = 0. 25, and thereafter f^ de- 
clines as increases. For trajectories with intermediate thrust, the f y function be- 
haves in qualitatively the same way. The peak value of f , however, decreases and oc- 

2 v 

curs at higher values of v w as ap Q is increased. For instance, with ap Q = 0. 01, 

V max = ^ a * v oo = 0*65. As a^ approaches the high-thrust regime, the peak be- 
comes low, very broad, and occurs at high values of v^. Although not shown on the fig- 
ure, fy — 1. 0 for all values of a^ when v^ — -1 (i. e. , when AVj m p — 0). 

The characteristic central angle 0 ^ is also influenced by v^ as illustrated in fig- 
ure 2(b) (p. 21). For impulsive thrust (ap Q = °°), 0^ is given by equation (16); it varies 

downward from a maximum of 180° at v^ = 0, asymtotically approaching a lower bound 
of 90° for large v^. For finite values of ap Q , the behavior is similar; the greatest 
value of 0 ( ,j 1 corresponds to v^> = 0, with a rapid approach to a lower bound as v^ in- 
creases. The larger values of 0 ^ at low v^, are due to the fact that low energy tra- 
jectories remain within the sphere of influence a longer time and thus are more strongly 
influenced by gravity. 

Effects of orbit acceleration parameter ap Q . - Although 0^ is a bounded function 

of v^ for given a po , it is not a bounded function of ap Q for given v w . This may be 
inferred from figure 2(b) and is further illustrated in figure 8, in which 0 ^ is plotted 
directly against ap Q for several typical maneuvers. The large values of 0^ at low 
a po are < * ue larg e value of 0^ Q associated with low-thrust spiraling. 

The behavior of fy as a function of ap Q is illustrated in figure 9 for several typical 
maneuvers. It is clear that f is bounded between the impulsive (f = 1) and ultra low- 

v v o 

thrust limits. The asymptotic value in the latter case is given in terms of v^ by equa- 
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Characteristic velocity ratio, f v - AV rh /AV 



Figure 8. - Effect of orbit acceleration parameter on trajectory deflection. Circular 
parking orbit; hyperbolic velocity parameter, 0.25. 



Figure 9. - Influence of acceleration and jet velocity parameters on gravity loss; hyperbolic velocity parameter, 0. 
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tion (25). It may be observed that f varies most rapidly with a in the r ang e 
0. 03 < a < 0. 3. This range includes the thrust to weight ratios that are normally of 

pu 

interest for vehicles using nuclear rocket propulsion and low parking orbits. 

The characteristic behavior of f may be understood by considering how the mean 

v 2 

tangential velocity component (v cos u) varies with a^ and v^. (For tangential thrust, 

only ( v) need be considered.) When equation (32) is used, the maximum value of f , 

2 ^ 
with given v^, occurs for minimum ( v) . Figure 10 depicts the velocity history of a 

typical set of tangential-thrust maneuvers (escapes from a circular orbit with Vj = 1) . 
The path velocity v during the thrusting phase is plotted against the cumulative mass 
ratio M./M bo on the logarithmic scale, and against Av cll on the linear scale. Each 
solid curve represents maneuvers with a given value of a^. The upper straight line 
represents impulsive maneuvers, while the lower curve denotes the zero-thrust limit. 
Loci of constant v w are shown in dashed lines. For high but finite acceleration, v in- 
creases monotonically with Av^, although not as rapidly as in the impulsive case. At 

low acceleration, however, v actually decreases initially, and only begins to increase 

o 

again after a fairly high value of v M has been attained. It is clear that for any final 



0 .5 1.0 1.5 2.0 

Characteristic velocity increment, Av ch , dimensionless 

Figure 10. - Velocity to mass ratio history for dimensionless escape. Jet velocity parameter, 1.0. 
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value of v^, <v) decreases as a^ decreases and this accounts for the fact that f y in- 
creases as a„„ decreases. 

po 2 

The location and magnitude of the peak value of f y as a function of v M may be pre- 
dicted by determining (v) (for example, by graphical means) from figure 10. For in- 
stance, in the low-thrust limit (lowest curve in fig. 10), it is readily verified that the min- 
imum value of (v) («0. 417)* occurs when Av ch « 1. 5 and v^ « 0. 5. These values in 
equation (4) yield 


f = 3 
v 


in agreement with both equations (25) and (32). 

o 

Effects of jet velocity parameter v-. - Variations in v^ and a have similar ef- 

J 


fects on escape and capture maneuvers. The effect of Vj, however, is qualitatively dif- 
ferent for the two types of maneuvers when they are compared on the basis of equal values 
of ap Q , as shown in figure 11. The reason for this behavior can be found by inspection of 
equations (9a) and (9b). That is, for finite Vj, a(T) increases continuously during an es- 
cape maneuver; the lower the values of Vj the greater the increase of a(r). As Vj — 0 
in an escape maneuver, f y therefore approaches the impulsive limit. On the other hand, 
a(r) decreases steadily during the backward integration of a capture maneuver so that as 
Vj — 0, f y approaches its low-thrust limit. For high values of Vj, a(r) for both types of 
maneuvers is nearly constant, and the values of f y approach a common asymptote. 

From a physical viewpoint, however, the difference between escape and capture ma- 
neuvers is more apparent than real. When escapes and captures are compared on the 
basis of equal initial acceleration, the values of f^ are similar. This may be recalled 
from figures 5 to 7 and is further illustrated by the lower dashed and solid curves in fig- 
ure 11. 


* Recall that (v) is found by averaging v with respect to Av^. Thus, in figure 10 

area under curve a - x 


< v >x = 


Av 


ch, x 


In the low thrust limit, the v - Av^ characteristic consists of straight line segments; 
then with v = v = 0. 5, for instance, 

X j X 

, v x _ 0. 5(1) (1) + 0. 5(0. 5)(1. 50 - 1) _ 0 , 

' v ' low thrust u - 41 ' 

Similarly, 


(v) = 1 

N / imp 


Av. 

+ .25 

2 
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Elliptic Orbits 


Elliptic planetocentric parking orbits are of interest because they offer a significant 
reduction in propulsive AV for a given value of v w , as suggested by letting e^ — 1 in 
equation (15). The reduction can approach 41. 4 percent of V c ref ; it is thus most signif- 
icant for low periapse orbits about the major planets. At Jupiter, for instance, the sav- 
ing can be as high as 10 miles per second, while at Mars it cannot exceed 0. 8 mile per 
second. 

Characteristic velocity ratios (f^) and geometric data (Vp Q and 0^) were presented 

in figure 3 for optimum-angle trajectories using an elliptic orbit with e_ = 0. 9. This 

po 

value of e po was selected as a compromise between AV savings (cf eq. (15)) and or- 
bital period considerations. 

It may be seen in figure 3 that the f data for elliptic orbits are similar in appear- 

* 2 

ance to those for circular orbits. The effects of v^, ap Q , and Vj follow the same qual- 
itative trends as in the circular orbit case. The two additional variables v and e 

po po 

introduced by the consideration of elliptic orbits are now discussed. 

Effects of initial true anomaly v . - For the optimum-angle trajectories considered 
P° 

herein, the power-on or -off point on the ellipse is simply chosen in order to mini- 

mize fy for given values of v^, ap Q , and Vj. The characteristic angle 0 ch is re- 
garded as a dependent variable; optimum values of 0^ (i. e. , corresponding to minimum 
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True anomaly on parking orbit, deg 


Figure 12. - Effect of true anomaly at elliptic parking orbit on escape 
or capture maneuver. Parking orbit eccentricity, 0.9; hyperbolic 
velocity parameter, 0. 25 ; jet velocity parameter, ®. 

f y ) are presented. The process of selecting i; to minimize f y is shown graphically 
in figure 12. The quantities f and 6 h are plotted against v for maneuvers where 

= 0. 9, v^, = 0. 25, and Vj = Impulsive and finite-thrust (a^ = 0. 1) maneuvers are 
represented by the solid and dashed curves, respectively. 

For the impulsive maneuver, the minimum value of f^ = 1. 0 occurs at the periapse 
(fp 0 = 0); this choice obviously yields the maximum value of (v) . The corresponding op- 
timum value of 0^ is 142°. The increase in f y for nonoptimum angles represents the 
effect of a geometrical constraint alone on the propulsive efficiency and mean velocity. 

For the maneuver with a^ = 0. 1, the minimum value of f y = 1 . 03 occurs at point a 
where i/ pQ = -47°. This choice maximizes (v) for equation (32) during the powered ma- 
neuver; it distributes the thrusting arc into roughly equal central angles ahead of and be- 
hind the original periapse. The associated optimum value of 9 ^ (point b) is about 143°, 
only slightly different from the impulsive value (142°). It may be noted that 0 ^ can be 
varied between 132° and 150° in this example for < 1, 10. 

Effect of parking orbit eccentricity e . - While v was selected in every case to 

po po 

yield minimum f v , the choice of e pQ is based primarily on mission objectives. Prob- 
ably the highest value of e^ would be used which is compatible with the observations, 
experiments, and operations that constitute the purpose of the mission. Although the 
value of e po = 0. 9 is believed to be a good compromise, it is quite possible that lower 
or higher values may be required in particular cases. 

The effect of other values of e__ is illustrated in figure 13 for the same maneuver 
2 po 

(v^ = 0. 25, a po = 0. 1, Vj = <*>) as was discussed previously. The minimum values of f y 
and the corresponding optimum values of ^ p0 and 0 ch are plotted against e . The 
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Figure 13. - Effect of parking orbit eccentricity on optimum tan- 
gential thrust maneuvers. Hyperbolic velocity parameter, 

0. 25; jet velocity parameter, dimensionless acceleration at 
parking orbit, 0. 1- optimum orbital true anomaly chosen for 
minimum characteristic velocity ratio. 


points a and b at e = 0. 9 represent the values of ^ and 0^ corresponding to the 
optimum i/ as shown in figure 12. 

It may be noted first that 0^ varies little with e^ Q . More important is the fact 
that f v decreases almost linearly with increasing e . Further examples (not illus- 
trated here) indicate that this behavior is typical for maneuvers where a^ > 0. 01. 
Therefore, in this medium- to high-thrust regime, values of f^ for eccentricities be- 
tween 0. 0 and 0. 9 can be estimated without major error from a linear formula such as 


v, e 


po 


= (1 


*> f v, 0. 0 + 


Xf v,0.9 


where 


x =^22 

0.9 


Comparison of Elliptic and Circular Orbit Results 

As was mentioned previously, the elliptic orbit data presented in figure 3 have the 
same general character as the circular orbit data of figure 2. There are, however, sig- 
nificant numerical differences. It was pointed out before that in the medium- to high- 
thrust regime (ap Q > 0. 01) f y tends to decrease with increasing ep Q . The reason for 
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this may be understood by recalling that the optimization of v amounted to placing as 
much of the maneuver as possible in a region of high velocity, that is, near the periapse. 
Hence, for short burn times (i. e. , high accelerations and/or low velocity increments), 

( v) ac j. « (v) j m p. For elliptic orbits, the velocity increment and hence the burn time are 
lower than they are for a circle; therefore, a lower acceleration can be tolerated before 
fy begins to increase sharply. 

Very low-thrust systems, on the other hand, require such long burning times that 
only a relatively small part of the maneuver can be placed near the periapse. For the 
first few revolutions at least, the greater proportion of time will be spent in the vicinity 
of the apoapse, and this results in a low value of (v). Therefore, as predicted by equa- 
tion (29) and shown in figure 3(a) (p. 30), the low-thrust values of f can become quite 
large as ep Q — 1. These high values cannot be reduced significantly by the use of vari- 
ational steering control. They are therefore an inherent feature of the trajectories being 
considered and cannot be attributed to the assumed tangential steering law. 

For these reasons it may be concluded that the AV savings predicted for elliptic 

parking orbits are actually enhanced by the finite-thrust effects if a„_. > 0. 01, but they 

po 

are degraded for very low-thrust systems. That is, the elliptic orbits are relatively 
more advantageous for high-thrust systems. 


APPLICATION TO MISSION PROBLEMS 

The foregoing discussion indicated the significance, interpretation, and use of the 
characteristic velocity correction factors (f y ) and geometric data that have been pre- 
sented. This section will illustrate, by typical examples, the application of the data and 
procedures already described to actual mission problems. 


Hohmann Transfers to All Planets 

The minimum AV transfer between planets which lie in circular coplanar heliocen- 
tric orbits is the well-known Hohmann (180°, contangential) trajectory. In view of equa- 
tion (3), the Hohmann transfers will also yield low propellant fractions. 

Table HI lists the pertinent trajectory parameters for Hohmann transfers from the 
Earth to the other planets (except Pluto) of the solar system. The planets are assumed 
to lie in circular coplanar heliocentric orbits, and the transfers begin and end in circular 
planetocentric parking orbits at 1. 1 planet radii. Representative values of I (900 sec) 
and Aj (0. 2 G 0 ) were assumed for the purpose of computing f y . The circular parking 
orbits at 1. 1 planet radii were chosen quite arbitrarily for the sake of comparison. 
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TABLE m. - IMPULSIVE VELOCITY CORRECTION FACTORS AND OTHER PROPERTIES 


FOR HOHMANN TRANSFERS FROM EARTH TO OTHER PLANETS 
[Specific impulse, I, 900 sec; initial thrust to weight ratio, A^, 0.2 G 0 .] 


Parameter 



Destination planet 




Mercury 

Venus 

Mars 

Jupiter 

Saturn 

Uranus 

Neptune 

- 

$ 

Dej 

9 

>arture £ 

a 

rom Earl 

% 

b 

6 

¥ 

Hyperbolic excess velocity 
at Earth departure, 

V M miles/sec 

4.66 

1. 56 

1.85 

5. 46 

6.38 

7.02 

7.26 

Circular velocity at Earth 
parking orbit, 

R po = *• 1R ® 

4.69 

4.69 

4.69 

4. 69 

4.69 

4.69 

4.69 

Gravity at Earth parking 
orbit, ft/sec 2 

29.6 

29.6 

29.6 

29.6 

29.6 

29.6 

29.6 

Impulsive velocity incre- 
ment to depart Earth, 
AV 0 , miles/sec 

3.42 

2.13 

2.19 

3.90 

4. 50 

5.00 

5.14 

Hyperbolic velocity param- 
eter, v^ 

0. 987 

0.111 

0. 155 

1.351 

1. 850 

2. 239 

2.391 

Dimensionless acceleration 
at parking orbit, a^ 

0.242 

0.242 

0.242 

0. 242 

0. 242 

0.242 

0.242 

Jet velocity parameter, Vj 
Characteristic velocity 
ratio, f y 

1.169 

1.169 

1.169 

1.169 

1.169 

1.169 

1.169 

1.107 

1.052 
Arrival i 

1.054 
it planet 

1.119 

1.132 

1.140 

1.142 

Travel time (Earth to 
planet), days 

106 

146 

259 

1000 

2200 

5900 

11 200 

Hyperbolic excess velocity 
at planet arrival, 

V c°,pl> mlles / sec 

5. 94 

1. 58 

1.65 

3. 51 

3.37 

2. 92 

2. 51 

Circular velocity at planet 
parking orbit, 

R po = 1 - 1R pl 

1.75 

4.29 

2.05 

25.2 

15.1 

8.89 

9.77 

Gravity at planet parking 
orbit, ft/sec 2 

9.6 

22.9 

10.65 

70.2 

. 31.2 

24.2 

27.4 

Impulsive velocity incre- 
ment to arrive at planet, 
AVpp miles/sec 

4.69 

1.94 

1. 29 

0.45 

6. 40 

3. 95 

4.28 

Hyperbolic velocity param- 
eter, 

11. 50 

0.368 

0.647 

0. 0194 

0.0498 

0. 108 

0.066 

Dimensionless acceleration 
at parking orbit, a^ 

1. 585 

0.406 

0.769 

0.69 

0.706 

0. 562 

0. 532 

Jet velocity parameter, Vj 
Characteristic velocity 
ratio, f y 

3.130 

1.277 

2.673 

0.218 

0.363 

0.616 

0. 561 

<1.01 

1.040 

1.021 

1.05 

(extrap- 

olated) 

1.052 

1.041 

1.042 



The hyperbolic excess velocity, the circular velocity, the local gravity, the impul- 
sive velocity increments, the dimensionless parameters, and f are listed for both ter- 
minals of each Hohmann transfer. For Earth, Mars, and Venus, f was found from fig- 
ures 5 to 7. For the capture maneuvers at the other planets, the previously described 

O 

iterative procedure based on equation (10) was used with figure 2 to find a^, v^, Vj, and 
thence f y . 

It is of some interest to note that with the chosen value of A .. (0. 2 G 0 ) the correction 
factors at the destination planets are only slightly greater than unity (e. g. , 1. 05 at Jupi- 
ter). In fact, they are smaller than for the associated Earth escape maneuver. This 
somewhat surprising observation is the result of two effects which together determine 

a po* 

(a) For capture maneuvers, a is related to a. by 

po 1 


a = 
po 


a i 


1 - m. 


= a.e 


AV ch/ G ® 1 


( AV ch/°® 1 ) 

For the maneuvers under consideration, the mass ratio e ranges from about 

2 (at Mars) to about 8 (at Jupiter). Thus A^ 0 , the acceleration at burnout, ranges from 
about 0. 4 to 1. 6 since A. was given as 0. 2 G^. 

(b) The ap Q is then related to A bo by equation (5c); that is, 


po 


^o 


G 


ref 


The gravity ratio G^/G re j ranges from about 1 to 3 for all the planets except Jupiter, 

where it is about 0. 4. For all planets, the net result is a fairly high value of a which, 

2 P° 

together with the generally low value of v^, leads to a low f y . 

Thus, for a one-way Hohmann trip to Jupiter, the neglect of finite-thrust effects 

leads to an error in AV ch of 0. 99 mile per second out of 14. 35 miles per second, or 

about 7 percent. 


Choosing the Circular Orbit Radius 

The low circular parking orbits previously assumed were sufficient for a compar- 
ative discussion of the behavior of at different planets. In an actual mission study, 
however, the choice of the parking orbit and the associated approach and departure ma- 
neuvers would receive very close attention. The low circular orbits have some definite 
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advantages, but it is also possible that more sophisticated orbits and maneuvers will be 
used in many cases. The purpose of this section, then, is to illustrate the use of the 
present data in finding an optimum parking orbit for a finite-thrust system. 

The simplest approach to parking orbit optimization consists of just selecting the cir- 
cular parking orbit radius Rp Q ^ that will yield the lowest value of AV^ for a given 
V^. Since the orbit is circular, the necessary hyperbolic direction for both the incoming 
and outgoing maneuvers (fig. 1(c), p. 6) can be attained with no AV penalty. Thus, the 
geometric effects are not considered in this example. 

Inspection of equation (6) reveals that AV- m p is a function of the parking orbit ra- 
dius alone for a circular orbit and given V^. It is readily seen that the minimum value 
of AVj m p occurs when 



or 


R 


2p 


po, opt 


pl 


(42) 


This is a well-known result for impulsive maneuvers. In the present case, however, it 
must be remembered that changing Rp Q , by changing V c ref and G ref , will affect the 
values of v^, a , and v.. Thus, the changes of f as well as of AV. must be 

J V JL111 JJ 

accounted for in seeking a minimum value of AV v* 

As an example, consider the problem of selecting the optimum parking orbit radius 


for a representative Earth/Venus round trip, which has the characteristics listed in 

table IV. This trip has a total duration of 450 days of which 30 are to be spent in the 

Venus-centered parking orbit. Data were ob- 
TABLE IV. - CHARACTERISTICS OF TYPICAL , . 

tamed from reference 1. 


450-DAY ROUND TRIP TO VENUS (1980) 


Maneuver 

Date 

(Julian) 

Hyperbolic excess 
velocity, 
EMOS 
(miles/sec) 

Earth departure 

2444330 

0. 131 (2.43) 

Venus arrival 

2444450 

. 165 (3.05) 

Venus departure 

2444480 

. 164 (3. 04) 

Earth arrival 

2444780 

. 298 (Atmospheric 



braking) 


Since the Venus arrival and departure ma- 
neuvers involve about the same value of V M , it 
is sufficient to consider only one of them, for 
example, the departure for which 
V^, = 3.04 miles per second. Figure 14 shows 
fy and AV ( ,| 1 as functions of R^/R^ and A. 
for this maneuver. 

In the lower portion, f is shown for values 
of R po /R^ ranging from 1. 1 to 16, I = 900 sec- 
onds, and initial thrust to weight ratios Aj of 
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°°, 0. 5, 0. 1, 0. 05, and 0. 01 G 0 . 

It may be noted that f falls off 
very rapidly as R increases; 

po 

even for A^ = 0. 05 G 0 , f y is incon- 
sequential when exceeds a few 

Venus radii. This, of course, is 
due to the fact that the term G ref 
in the equation (5c) defining 

a po (= A po /G ref ) falls off as 1/R po’* 
hence for a given value of A. , a^ 

i pu 

will become large at a sufficiently 
high radius. 

The upper curve of figure 14 

shows the values of 

AV_v, = f AV. that result when 
ch v imp 

the effects of R p0 on both f y 
are taken into account 


Figure 14. - Optimization of circular parking orbit radius for Venus 
escape maneuver. Hyperbolic velocity, 3. 04 miles per second; 
specific impulse, 900 seconds. 


and 

AV-_^ are taken into account. The 
imp 

lower solid curve represents impul- 
sive maneuvers (Aj = °°); as pre- 
dicted by equations (33) and (34), the 

minimum value of AV. mp (2. 16 miles/sec) occurs where Rpg/Rg = 4. 4. This radius is 
also nearly optimum for values of A i as low as 0. 1 G 0 ; for smaller values, however, 
the optimum radius is higher. For instance if A^ = 0. 01 G 0 , then the minimum value of 
AV ch (2.25 miles/sec) occurs at a radius R^ ~ 10 R^. This is a lower characteristic 
velocity than the low orbit (R^/Ro = 1.1) impulsive value (2. 48 miles/sec), yet the low 

po If 

value of Aj implies that a very small and light propulsion system could be used. 

This example not only illustrates the effect of finite thrust on parking orbit selection, 
it also demonstrates a possible area of application for propulsion systems (such as arc 
jets) which might otherwise appear unattractive. 


Optimization of Initial Thrust to Weight Ratio 

Another frequently occurring mission problem that can be solved with the aid of the 
f y data presented here is the selection of the optimum value of A^. This is an important 
area of study, not only from the standpoint of minimizing the vehicle gross weight, but 
also because it has a direct bearing on the problem of selecting engine sizes for future 
development. For given values of V w and R pQ , a high value of A^ implies a low value 
of m p but a large, heavy engine: a low value of A i implies a light engine but a high 
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Figure 15. - Vehicle weights for typical Earth escape maneuver. 

Hyperbolic velocity, 2.43 miles per second; parking orbit radius, 

1.1; specific impulse, 900 seconds; propellant sensitive structure 
fraction, 0.1; specific weight of engine and thrust structure, 

0.2; acceleration sensitive structure fraction, 0.01. 

propellant fraction. Thus, it may be anticipated that minimum M will occur for inter- 
mediate values of A^, when there is a proper balance between propellant- and engine- 
related weights. 

This process is illustrated by using the Earth-departure maneuver in table IV as a 
final example. In figure 15, the initial space vehicle gross mass in a low circular Earth 
parking orbit is plotted as a function of the initial acceleration. The propellant fractions 
are computed by using equation (3) and figure 5(a) (p. 42), with I = 900 seconds. Gross 
weights are then estimated from 


M - M. 1 * A max k as 

(1 - - m p k ps - A,k ts 

where the following arbitrarily assumed parameter values were used: 
Total Earth-departure payload, including upper stages: 

= 500 000 lb 

Propellant sensitive structure fraction: 


(43) 


M 

kps =— E£=0.10 
P 8 M P 
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Specific weight of engine and thrust structure: 



0. 20 


Acceleration sensitive structure fraction: 


M. 


k as = 


as 


A max M f 


0.01 


The upper curve representing has a minimum value of about 900 000 pounds when 
= 0. 225 G^. The corresponding nuclear engine would have a thrust rating of about 
200 000 pounds. The lower curves, which denote, respectively, the propellant and tank- 
age weight, engine and thrust structure weight, and interstage structure weight, indicate 
how M is distributed among these major weight items. 

o 


CONCLUDING REMARKS 

A wide range of data defining the essential dynamic and geometric characteristics 
(fy and 0 ^) of continuous finite-thrust trajectories has been presented. These results 
are presented in dimensionless form and may be easily scaled to represent any case of 
physical interest. In this way, accurate propellant fractions and realistic trajectory ge- 
ometry may be conveniently derived from impulsive interplanetary trajectory data. 

The characteristic velocity ratio fy is interpreted as a measure of the mean propul- 
sive efficiency ri^ r m along the powered arc. Both f v and 77 ^ m are shown to be ex- 
plicit functions of the mean tangential path velocity component (v cos u) . This quantity, 
and hence also f y and ri^ r m , obviously are directly affected by the steering program 
u(t) that is used, whether a significant gravity field is present or not. Thus, in the gen- 
eral case, f y is as much due to steering losses as to gravity losses. 

Of particular interest in this report is the class of optimum-angle trajectories. 

These are so named because both u(r) and the initial power-on point are chosen only to 
minimize fy (by maximizing (v cos u)) without regard to the consequent value of 0 ^. 

For such trajectories, it is shown in the appendix that tangential steering yields very 
nearly minimum values of fy, especially in the case of an elliptic orbit. 

These optimum- angle trajectories can always be used with circular parking orbits, 
because any final direction of can be attained without penalty in f by merely selec- 
ting the appropriate power-on point. Optimum- angle trajectories can also be used for any 
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one maneuver based on an elliptic orbit by orienting the major axis at an angle 0^ to 
V*. When two maneuvers are based on the ellipse, however (that is, during the capture - 
orbit- escape sequence at the destination planet of a round trip), it is generally necessary 
to reorient the ellipse (by means of auxiliary maneuvers) before optimum-angle trajecto- 
ries can be used for both escape and capture. It remains to be seen whether or when the 
steering losses due to optimum constrained angle trajectories would result in a smaller 
overall penalty than the auxiliary orbit shifting maneuvers. 

By its definition f y is bounded above 1 . An upper bound is also derived in the very 
low-thrust limit. This upper bound increases monotonically as e^ increases, ranging 
from 3. 0 for a circular orbit to 10. 6 for e = 0. 9, and becomes very large as e —1.0. 

rs* PO 

Thus for very low-thrust systems (a < 0. 001), f y is considerably larger for a maneu- 
ver from an elliptic orbit than it is for an equivalent circular orbit maneuver. This par- 
tially offsets the AV advantage of elliptic orbits. That is, the saving in the low-thrust 
AV^ is smaller than the savings in AVj m p. 

In direct contrast, medium- to high-thrust systems (ap Q % 0. 01) have lower values of 
fy for elliptic orbits than they do for equivalent maneuvers from a circular orbit. There- 
fore the advantage due to elliptic orbits is actually enhanced in this case (in the sense that 
the saving in characteristic velocity for finite thrust is greater than the saving in AV^ m p). 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, April 12, 1966, 
121-30-02-01-22. 
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APPENDIX - COMPARISON OF TANGENTIAL AND OPTIMUM 
STEERING CONTROL FOR ESCAPE FROM CIRCULAR 
AND ELLIPTIC PARKING ORBITS 

Herein the tangential and optimal steering controls are compared for the problem of 
constant- thrust escape from circular and elliptic orbits. Previous work (reported in 
ref. 9) has shown that in the absence of geometric constraints, tangential steering is 
nearly optimal for constant acceleration departure maneuvers from a circular orbit to 
parabolic final energy (v^ = 0) . These results are now extended to include maneuvers with 
v oo ^ constant thrust, and elliptic parking orbits. 


The Maximum Principle 

The optimal angle- of- attack control program u(r) is defined by Pontryagin f s maxi- 
mum principle (described in ref. 10). According to this, the control u(r) can be optimal 
only if the function 


= i^r' + ^ 2 V ' + + ^5 Av^ 


(Al) 


attains its maximum with respect to u(t) for 0 < r < t^ q . Furthermore, M is a constant, 
and its value is zero for an optimum-time trajectory. The state variable derivatives are 
given by equations (7) and (8), and the adjoint variables are defined by 


1 3r 


2 sin a . ( 2 v\, /v \ , 

^ 2 3 7 cos 7 4: 


sin of +i//. 

a po sin u /l 1 \ 

- — + (- + icos a 

1 

l-m_ 2 \r 2 2/ 

L p v \ r v / 


cos a 


=-Wco» a+ *,M2L»W J z - 1 


da 


— | (- sin a) + x[/ A — (- sin a)\ 
rV 4 r I 


V r r 2 v 


y (A2) 


&e 


ip], =0 or ip, = constant 

* dd q 


n = - 


M 


d Av 


= 0 or i|/^= constant 


ch 
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Since the range of u(t) is not restricted, the maximum condition reduces to 




i^ = o 

3u 


au 2 


< o 


and this leads to the optimal steering control law 


tan u(t) = 


\Jj 0 v 


(A3) 


(A4) 


Transversality Conditions 


Ten boundary conditions are required to specify the simultaneous solution of equa- 
tions (A2) together with the equations of motion (7) and (8). For initial conditions are 
specified in terms of by equation (12), and one terminal condition is given in equa- 
tion (13). Since the xf/'s are homogeneous first-order functions, one of the initial values 
can be taken as a scale factor. The remaining four relations are defined by applying the 
transversality condition (ref. 10) at both the initial and final times. 

At the initial time r = 0, 


sin = 

P° . 2 


cos^ a 


' cos v + e pol 


\p t r +^ 2 /v L 


1 + e. 


po 


1 + e po^4 
6 po J 


(A 5) 


This condition states that if a is normal to the one-dimensional manifold defined by equa- 
tion (12). 

At the terminal time, the transversality condition yields the final three relations: 


^3, bo “ ^4, bo " 0 


(A6) 


^2, bo 

^l,bo 


= r 2 V 
r bo v bo 


(A7) 
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Solutions 


The two-point boundary value nature of the problem requires an iterative solution to 
search in values of ^(0) that will lead, by means of equations of motion (7), the adjoint 
equations (A2), and the maximum condition (A4), to the desired terminal values (13), (A 6 ), 
and (A7). Although five initial values of 1 ^( 0 ) are needed, all but two of these can be ob- 
tained without a numerical search. 

(1) After j // 1 . - ip. . have been selected, (a constant) is chosen so as to make 

1, X ft, 1 o 

<# = o. 

(2) In view of (A 6 ) and (A2), 1^4 = 0. 

(3) Then, the initial value of is set equal to unity as the scale factor. 

Thus, it only remains to find i ^3 j. For elliptic orbits, these may be found in 
terms of assumed values of u. and by inverting equations (A4) and (A5) : 



Finally, the initial values u. and v are used as the independent variables in a two- 
dimensional search to satisfy the final conditions (13), (A 6 ), and (A7). The use of these 
bounded, physically significant variables in place of the abstract ip's results in a consid- 
erable practical simplification of the problem. The values of u.(0) and u that corre- 
spond to an optimum tangential trajectory (figs. 3(s), (t), (u), and (v) and 4 to 8 ) provide 
excellent starting points for the search. 

Conditions (A 6 ) and (A7) have the further implication that for optimum-angle trajec- 
tories (where ^4 = 0) both ^3 and its first derivative ^3 vanish at Thus, 

in view of equation (A4), the final portion of an optimum-angle trajectory is tangential to 
a high order. Any improvement due to optimal steering must therefore be generated in 
the early part of a trajectory. In the high-thrust limit, ^ remains bounded and there- 
fore 1^3 is continuous across an impulse. That is, tangential steering is the optimum 
control for impulsive thrust. Further, the transversality condition (A5) shows that this 
impulse is to be applied at periapse (y = 0 ) since . = 0 . 
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TABLE V. - COMPARISON OF TANGENTIAL AND Numerical ReSUltS 

It is clear that tangential steering 
cannot satisfy the conditions of the 
maximum principle when finite thrust 
and nonzero burning times are consid- 
ered. Nevertheless, there are two 
reasons to expect that tangential steer- 
ing would yield nearly minimum Av^: 

(1) As pointed out before, at least 
the latter part of an optimal trajectory 
must be tangential to a high order. 

(2) Tangential steering, as is well 
known, yields the maximum rate of in- 
crease of at any instant. In 

table V, optimal- and tangential- thrust maneuvers are compared for escape from a cir- 
cular orbit with ap Q = 0. 1. It is clear that in no case does optimal steering reduce f y 
by more than 1 percent below the tangential value. Thus, the result of reference 9 is ex- 
tended to the case of constant thrust as well as constant acceleration and to values of v 

OO 

in excess of zero (parabolic energy). 

The effect of the eccentricity of an elliptic parking orbit is next illustrated in ta- 
ble VI. Here it is evident that, for the conditions chosen, the difference in f v between 
tangential and optimal steering is greatest for circular orbits and decreases rapidly as 


OPTIMAL ESCAPES FROM CIRCULAR ORBIT 
[Dimensionless acceleration at parking orbit, a^, 0. 1. ] 


Hyperbolic 

Jet velocity 

Characteristic velocity ratio, f y 

velocity 

parameter, 

V* 

parameter, 

v i 

Tangential escape Optimal escape 

0 

OO 

1. 250 a l. 241 

0 

1.0 

1. 179 1. 171 

.10 

1.0 

1. 199 1. 190 

.25 

1.0 

1.225 1.214 

. 50 

1.0 

1. 261 1.250 

1.00 

1.0 

1. 300 1. 290 


a From reference 9: all other values from present analysis. 


TABLE VI. - EFFECT OF PARKING ORBIT ECCENTRICITY 

ON CH ARACT ERISTIC VELOCITY RATIO FOR 

TANGENTIAL AND OPTIMAL STEERING 

o 

[Hyperbolic velocity parameter, v^, 0. 1; dimensionless 

acceleration at parking orbit, a„_, 0. 1; jet velocity 

po 

parameter, °°. ] 


Parking orbit 

Characteristic velocity ratio, f y 

eccentricity, 

e po 

Tangential steering Optimal steering 

0. 9 

1 . 011 1 . 011 

.8 

1. 024 1. 024 

.6 

1.066 1.065 

.333 

1. 150 1. 147 

.00 

1. 270 1. 260 
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TABLE VH. - EFFECT OF ORBIT ACCELERATION PARAMETER 


ON CHARACTERISTIC VELOCITY RATIO FOR TANGENTIAL 
AND OPTIMAL STEERING 


[Parking orbit eccentricity, e , 0. 9; hyperbolic velocity 
parameter, v^,, 0.1; jet velocity parameter, «. | 


Orbit acceleration 

Characteristic velocity ratio, f 

parameter, 

a po 



Tangential steering 

Optimal steering 

0. 1 

1.011 

1.011 

.01 

1.481 

1.480 

.001 

4. 261 

4.256 

.0001 

6. 890 

6. 871 


increases. That is, tangential steering is a better and better approximation to opti- 
mal as ep Q increases. 

Finally, in table VII, the effect of a.^ o on the comparison of tangential and optimum 
steering is illustrated for an elliptic orbit. From this it may be seen that, for all values 
of ap Q , the difference in f y for optimum as compared to tangential steering is a small 
fraction of 1 percent when an elliptic orbit is used. This is true even for very low accel- 
erations for which is large. Thus, the high values of f y predicted (e. g. , by 
eq. (29)) for elliptic orbit departures with low ap Q and v^ are a true and inherent fea- 
ture of the trajectories being considered; they are not due to any peculiarity of the tan- 
gential steering law. 

In summary, it may be inferred from the foregoing discussion that across-the- 
board incorporation of optimal steering control would not modify the conclusions nor sig- 
nificantly change the data presented in the body of this report. 
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